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Abstract 

We study the evolution of heavy quarkonium states with temperature in a Quark Gluon Plasma (QGP) 
by evaluating the in-medium Q-Q T— matrix within a reduced Bethe-Salpeter equation in both S— and 
P— wave channels. The underlying interaction kernel is extracted from recent finite-temperature QCD lattice 
calculations of the singlet free energy of a Q-Q pair. The bound states are found to gradually move above the 
Q-Q threshold after which they rapidly dissolve in the hot system. The T— matrix approach is particularly 
suited to investigate these mechanisms as it provides a unified treatment of bound and scattering states 
including threshold effects and the transition to the (perturbative) continuum. We apply the T— matrix 
to calculate Q-Q spectral functions as well as pertinent Euclidean-time correlation functions which are 
compared to results from lattice QCD. A detailed analysis reveals large sensitivities to the interplay of bound 
and scattering states, to temperature dependent threshold energies and to the "reconstructed" correlator 
used for normalization. We furthermore investigate the impact of finite-width effects on the single-quark 
propagators in the QGP as estimated from recent applications of heavy-quark rescattering to RHIC data. 

PACS numbers: 25.75.Dw, 12.38.Gc, 24.85. -l-p, 25.75.Nq 
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I. INTRODUCTION 



Bound states of heavy quarks (charm and bottom, Q = c, b) have long been recognized as 
valuable objects for spectroscopy in Quantum Chromodynamics (QCD), thereby illuminating the 
nature of the static quark-antiquark potential (cf. Ref. for a recent comprehensive overview). 
This opportunity carries over when embedding quarkonia into hot and/or dense matter, providing a 
rich laboratory for the study of medium modifications. The latter include (Debye-) color-screening 
of the Q-Q interaction, dissociation reactions induced by constituents of the medium, and the 
change in thresholds caused by mass (or width) modifications of open heavy-flavor states {D and 
B mesons or c and b quarks). The challenge is to develop a theoretical framework that allows a 
comprehensive description of heavy quarkonia in the Quark-Gluon Plasma and their production in 
ultrarelativistic heavy-ion collisions. 

Lattice QCD (IQCD) calculations have made substantial progress in characterizing in-medium 
quarkoniumproperties from first principles. In particular, it has been found that ground state char- 
monia [2, 0, |j] and bottomonia [5] do not dissolve until significantly above the critical temperature, 
Tc. This finding has been qualitatively supported in model calculations based on potentials ex- 
tracted from IQCD, using either a Schrodinger equation to solve the bound stateproblem 0,0,05 0]) 
or a T-matrix approach which simultaneously accounts for scattering states [lO[]. The survival of 
low- lying quarkonia above Tc, in connection with effects of color-screening, parton- induced disso- 
ciation and medium modified open-charm and -bottom thresholds, has recently been implemented 
for heavy-ion collisions [Til . 

A more quantitative (and reliable) comparison of model calculations to IQCD can be performed 
at the level of (spacelike) Euclidean-time correlation functions [lil . The latter are directly 
evaluated in IQCD with good accuracy, while the conversion of (timelike) spectral functions as 
evaluated in model approaches merely involves a straightforward convolution with a thermal weight 
function (as opposed to an inverse integral transform when going from Euclidean to Minkowski 
space). One of the challenges in such studies is that the model calculations need to describe 
not only the bound-state part of the spectral function but also its continuum part as well as 
threshold effects. In Ref. [9] a quantitative calculation of Euclidean correlators was performed using 
temperature-dependent heavy-quark potentials in a Schrodinger equation. The latter has been used 
to determine the bound-state spectrum in (5-function approximation (characterized by a binding 
energy and amplitude (or decay constant)), while the (onset of the) continuum was approximated 
with perturbation theory. While general trends of the IQCD correlators were captured, significant 
discrepancies were established especially in the 5- wave charmonium channels (r/c and J/ip). In 
particular, the importance of a reliable treatment of the continuum threshold was recognized. 

In the present paper we evaluate charmonium and bottomonium correlators using a different 
method. The basic input are still in-medium Q-Q potentials as estimated from IQCD, but we 



will employ these within a scattering equation to calculate the in- medium Q-Q T— matrix [10|]. 
The main advantage of the T— matrix approach is that it simultaneously incorporates bound and 
scattering states based on the same interaction. Especially for situations of dissolving bound states 
(as expected for the problem at hand) the T— matrix provides a more comprehensive, and thus 
more reliable, description of the underlying nonperturbative effects. At the correlator level, the 
high-energy limit can be recovered by appropriate normalization of the uncorrelated (perturbative) 
limit, and no decomposition into bound-state and continuum parts is necessary. In addition, 
the T— matrix equation allows for a straightforward implementation of in-medium single-particle 
(quark) properties via pertinent self-energy insertions in the two-particle Green's function, which 
we will also investigate. 

Our article is organized as follows: in Sec. |TI] we recall the basic set-up of, and input to, 
the two-body scattering equation, including partial-wave expanded potentials and single-quark 
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selfenergy insertions. In Sec. lIIII we evaluate the finite-temperature T— matrices for S— and P— wave 
quarkonia; we first extract heavy-quark potentials from IQCD in Sec. lIII Al including a discussion of 
its short- and large-distance limits and relations to single-quark properties, followed by our baseline 
results for the finite-temperature quarkonium T— matrices in Sec. IIIIBI In Sec. IIVI the latter are 
first employed to construct pertinent spectral functions (Sec. IIV Ap . followed by a calculation of 
Euclidean correlators (Sec. IIVBI) and a discussion of their properties in comparison to other model 
and IQCD results (Sec. IIV C]) . Sec. IVl contains our conclusions and an outlook. 



II. SCATTERING EQUATION AND IDENTIFICATION OF BOUND STATES 

We here summarize the main features of the T— matrix approach to study quark-antiquark in- 
teractions in the QGP, as employed in Ref. It utilizes a three-dimensional reduction of the 
Bethe-Salpeter equation which neglects virtual particle-antiparticle loops and amounts to resum- 
ming the scattering series in ladder approximation. The pertinent Lippmann-Schwinger equation 
for the off-shell T— matrix in a given partial- wave channel (specified by angular momentum /) 
reads^ 

TiiE;q',q) = Vi{q',q) + -J^ dk ^ Vi{q' , k) Gqq{E; k) Ti{E- k, q) [1 - 2/"^ {ut)] , (1) 

where q {q') are the incoming (outgoing) relative quark three-momenta in the center of mass 
(CM) frame and E is the CM energy^. Eq. ([T]) is written for vanishing total 3- momentum of the 
heavy-quark pair, which gives the above (simple) form of the Pauli blocking factor with f^{uj) = 
[exp(u;/r) -|- The intermediate two-particle propagator is evaluated in the Blankenbecler- 

Sugar reduction scheme [l^ (uncertainties due to other reduction schemes have been checked to 
be small [l(j|). 

TUp 1 

GgqiE; k) = — ^/A - col - 2iu;kluiJ:{u;k,k) ' 

where uj^ is the solution of the quark dispersion relation, 

ujk = V "rm? + k"^ + Re S(a;fc, k) , (3) 

with a quark-mass term (m) and selfenergy (S) to be discussed below. The interaction kernel of 
the scattering equation, Vi{q',q), is provided by the heavy-quark potential in momentum space. 
It follows from a Fourier transformation of the coordinate-space potential, V{r), which we obtain 
from IQCD calculations as elaborated in Sec. IIII Al below. The components of the potential in the 
partial-wave basis are given by 

Vlici.q) = ^ j^\u,.,V{q\q)Pi{u,,,) = i- £\u,,,Pi{ug,,) J d^rV{r)e'^'^-'^'^^^ , (4) 

with Pi{x) the Legendre polynomial of degree I and Uq/q = cos < q,q' >. 

The T— matrix equation ([1]) is solved with the algorithm of Haftel and Tabakin 16|]: after 
discretizing the momentum integration, Eq. ([T|) is converted into a matrix equation, 

N 

Y,T{E),kT{E)uj = V,j , (5) 
fe=i 



^ The partial wave expansion reads T = 47r + 1) Ti Pi{cos6), and similarly for the potential. 

^ The on-shell T— matrix is defined for q = q and E — -^s = 2ijjq, with uiq the (relativistic) on-shell heavy-quark 
energy. 
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where, schematically, !F = 1 — w V Gqq [1 — 2/*^] (with w denoting an integration weight). The 
solution for the T— matrix then follows from matrix inversion. 

To assess the presence of heavy quark-antiquark bound states, the T— matrix has to be studied 
below the Q-Q threshold, = 2a;g=o. The non-relativistic potential, Vi{q' , q), is only defined for 
real external three-momenta, and therefore an evaluation of T below the Q-Q threshold requires a 
prescription for the subthreshold continuation of the potential. For S"— wave scattering we follow 
the standard convention of setting the momenta to zero [l7| , 



To{E < Eth) = To{E; q' = q = 0) . (6) 

The reliability of this continuation can be checked by exploiting the (numerical) matrix form of the 
scattering equation. Since a bound state corresponds to a pole of the amplitude on the real energy 
axis below threshold^ follows that the determinant of the transition matrix T must vanish at the 
bound-state energy [la], 

detTiE) = 0, E<Eth. (7) 
A similar condition arises from the solution of the Schrodinger equation for the bound state prob- 



lem 16|, ll8|]. This is equivalent to finding the zeroes of the Jost function in scattering (S— matrix) 
theory. For P— wave states, the potential is proportional to the in- and outgoing quark momen- 
tum and therefore the continuation in Eq. ([6]) cannot be applied. However, the condition Eq. ([7]) 
remains valid and will be used to determine P— wave bound states. 

The quark selfenergy figuring into the two-particle propagator, Eq. receives contributions 
from interactions with (light) anti-/quarks and gluons in the heat bath. In Ref. flo| this was 
schematically written as 



j: = j: + J fTQgSg, (8) 

where S denotes the gluonic piece and the second term involves the heavy- light quark T— matrix 
closed by a light-quark pro pag ator, Sq, and a thermal distribution, f^. Rather than using an explicit 



model calculation for Tg^ jl9t | , in the present work we will constrain ourselves to the following levels 
of approximation: (a) a fixed heavy-quark mass m (i.e.. Re S = 0) together with a small imaginary 
part, ImS = —0.01 GeV, mostly for numerical purposes (to avoid 5-function like bound states in 
the T— matrix); (b) a temperature dependent heavy-quark mass as estimated from the asymptotic 
value of the IQCD heavy-quark internal energies; (c) a heavy-quark width as calculated in an 
effective model for resonance (plus perturbative gluon) interactions in the QGP [20], which has 
been shown to give reasonable agreement with data on suppression and elliptic flow of semileptonic 
electron spectra from heavy-quark decays in Au-Au collisions at RHIC [21|]. We also note that 
interactions with heavy anti- / quarks from the medium can be safely neglected due to the smallness 
of the number of Q's in the system. This is different to (and simpler than) the situation of the 
light-quark selfenergy which, in turn, figures into the calculation of the T— matrix, constituting a 
self-consistency problem as has been evaluated, e.g., in Ref. [Ifli]. 

III. TEMPERATURE EVOLUTION OF HEAVY QUARKONIUM T-MATRICES 
A. Quark- Antiquark Potential from Lattice QCD 

For the driving kernel of the scattering equation we focus on heavy-quark potentials from IQCD. 
There is an ongoing discussion as to how to properly extract them from the (static) heavy-quark 
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free energy. Using directly the latter leads to a dissociation temperature of ground-state charmonia 



of about T^iss — 1-1 Tc [24], while the lattice analysis of spectral functions suggests that rjc and 
J ftp survive up to around 2Tc. These higher values for Tdiss can be recovered if the (color-singlet) 
internal energy, 

C'.^F.-r^, (9) 

is identified with the Q-Q potential 0, B@, Another possibility, namely a suitable linear com- 
bination of Ui and Fi, has been suggested in Ref. The temperature derivative of discrete IQCD 
"data points" involved in the extraction of Ui induces significant uncertainty which is comparable 
to, e.g., the difference between quenched and unquenched results (after rescaling of the critical 
temperature), as studied in Ref. [10(]. In view of this situation, one has to accept a certain level of 
uncertainty in the potential. To adequately account for this, we adopt two versions of the internal 
energy, Eq. (|9]), as the potential: (i) based on fits to the two-flavor IQCD results for the free energy 



from Refs. [23l . |24| |. we explicitely perform the temperature derivative in Eq. Q; (ii) we directly 
fit the two-flavor IQCD internal energy data as extracted in the calculation of Ref. [25(. Further 
investigations of the impact of using different definitions of the Q-Q potential will be considered 
in future work. 

The long-distance limit of the potential entering Eq. ([T]) has to be normalized to zero to ensure 
the convergence of the scattering equation, 

V{r,T) = Ui{r,T)-U^{T) , (10) 

with C/f°(r) = Ui{r^ 00, T). In Refs. 0, B S 0, l^l the hnearity of the Schrodinger equation is 



exploited to trade the internal energy at infinite distance into the energy of the bound states. In 
Ref. B, C/f (T) is interpreted as an effective in- medium contribution to the quark mass, AmQ{T) = 
[/f°(T)/2 and implemented as a change in the Q-Q threshold energy, Efh = 2mQ + C/f°(T), in the 
calculation of the mesonic spectral function. It is argued that this correction should not modify 
the mass operator in the Schrodinger equation since quarks inside a bound state do not "sense" 
the medium and therefore should not be subject to medium-induced mass modifications. This 
ambiguity in the interpretation of the internal energy at infinite distance can be resolved within 
the many-body scattering equation approach. The interaction of the quark with the surrounding 
medium induces a selfenergy which is encoded in an effective mass change (Amg = ReE). As 
such, this medium effect has to be included in the two-particle propagator, Eq. ([2]), and therefore 
contributes in a nonlinear way as it is iterated to all orders in the scattering equation series. As 
we discuss in the following sections, this not only modifies the Q-Q threshold energy but the 
evolution of the binding energy (total mass of the bound state) with temperature, leading to 
different dissociation temperatures when the effective in-medium mass is considered. 

Another intriguing problem regarding the interpretation of C/f° as an in- medium quark mass (or 
selfenergy) is its possible momentum dependence. Whereas at infinite separation (low momentum 
transfer) operates as Amg, at short distances (high momentum transfer) the Q-Q system 
is no longer sensitive to the medium, and the quark mass should be unaffected. Implementing 
such a momentum dependent quark selfenergy would require a microscopic treatment of the mass 
subtraction of the internal energy. Instead, we will consider two limiting scenarios, namely (a) no 
in-medium mass correction and (b) in-medium effective mass as given by AmQ(r) = U^iT) /2. 

Finally, for the potentials in momentum space we introduce a relativistic correction motivated 



by the velocity- velocity (Breit) interaction in electrodynamics 28|], see also Ref. [10.] . It amounts 
to the following factor: 

V{q\ q) - y(g', q) [1 + g'V^,'']'/' [1 + ■ (H) 
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The two potentials used in this work are summarized in Fig. [TJ For case (i) discussed above 
(extraction from the IQCD free energy), it evolves rather smoothly with temperature (left panel), 
while for case (ii) (IQCD internal energy) the potential is initially more attractive but weakens 
rapidly with temperature (right panel) for T ^ 1.5 Tc and slows down thereafter. The rapid 
decrease is mostly induced by ?7f°(T), see also Ref. [25 1. 




B. Quarkonium T— Matrices in the QGP 

We now turn to the numerical results for the finite-temperature T— matrices in the c-c and 
h-b sectors, obtained by solving the scattering equation ([1]) in both S— and P— wave channels as 
described above. 



1. S-Wave States 



In a first step, we consider the case of narrow quark spectral functions with ImS = —10 MeV 
(for numerical purposes) and constant (temperature- independent) heavy-quark masses (ReS = 0). 
The latter are fixed so that the corresponding ground states are located approximately at their 
vacuum masses for the lowest considered temperature (T = 1.1 Tc), yielding rric = 1.7 GeV and 
rub = 5.15 GeV. 



Fig. [2] summarizes the on-shell S"— wave c-c scattering amplitude as a function of CM energy, 
for several temperatures from 1.1 Tc to 3.3 Tc, as well as the determinant function detJ-{E) (in 
arbitrary units). Since we do not include the hyperfine (spin-spin) interactions, r]c (rjb) and J/ijj 
(T) states are degenerate. At the lowest temperature, we recover the charmonium ground state 
at ~ 3.0 GeV, and also find a cusp at the cc threshold energy indicating that the first excited 
state (ascribed to the ip') has just melted. The determinant function, det^(-E), vanishes exactly 
at the ground state energy coinciding with ReT(£^) = 0, thus corroborating our subthreshold 
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FIG. 2: Real and imaginary parts of the T— matrix for 5— wave c-c scattering in the QGP based on potentials 
derived from the IQCD free energy of Ref . [24| . Also shown is the determinant function det J- (dashed line, 
arbitrary units). From left to right and up to down the temperatures are (1.1, 1.5, 2.0, 2.5, 3.0, 3.3) Tc- 



continuation of the T— matrix, Eq. For higher energies det^(£') approaches zero again, but 
as soon as the threshold is reached it deviates indicating that the first excited state has already 
crossed into the continuum spectrum^. As the temperature is increased, the bound charmonium 
state gradually moves toward threshold, indicating a reduction of its binding energy. At the same 
time, the magnitude of the T— matrix is appreciably reduced. The J/'ip{lS) survives as a bound 
state well beyond Tc, eventually crossing the threshold at about (2.8 — 3.0) Tc, after which it turns 
into a resonance and rapidly melts in the hot system^. Our results agree reasonably well with 
those of Ref. [iB], with some differences in size and shape of the scattering amplitude, in particular 
a larger dissociation temperature. This is mostly due to a different parametrization of the Q-Q 
potential (reflecting the uncertainties in a derivation of the Q-Q internal energy from a fit to the free 
energy, cf . Sec. HUB 41 below) and a different choice of the two-particle propagator, which implies 
deviations at order 0{p/m). The robustness of this approach to dynamically generate quarkonium 
bound and scattering states and their evolution with temperature is confirmed and extended in 
the following to study other quarkonium states. 



The results for 5— wave b-b scattering are depicted in Fig. [3l At the lowest temperature the 
T— matrix exhibits two bound states, as well as the remnant of a third one. The bound-state 
locations are again quantitatively confirmed by the vanishing determinant of the transition matrix 



^ Strictly speaking, J-{E) is not a purely real function since we have included a small imaginary part in the two- 
particle propagator. More precisely, the bound-state condition, Eq. reads Re{det jr(iJ)} = 0. We are plotting 
Ke{AetT{E)} in Figs. [1 13] and H 

* We refer to a state "melting" or "dissolving" when the scattering amplitude is strongly broadened and diminished 
corresponding to a loss of the resonant structure. 
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FIG. 3: Same as in Fig.[2]but for 5— wave b-b scattering. From left to right and up to down the temperatures 
are (1.1, 1.5, 1.8, 2.1, 2.7, 3.5) Tc. 
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FIG. 4: Pole identification function (det T) for charmonium (left) and bottomonium (right) P— wave scat- 
tering at several temperatures. 



(dashed lines), while it barely reaches zero at the location of the third structure in the T— matrix, 
which carries much smaller strength, indicating that it has practically melted in the medium. The 
two bound states at E ^ 9.35, 10.05 GeV are ascribed to the ground and first-excited bottomonium 
states T(15), rjb and T{2S), r/^, respectively. The T{2S) moves above the b-b threshold at T sa 
1.8 Tc, whereas the IS* state survives in the QGP until much higher temperatures, beyond T « 
3.5 Tc. 



2. P- Wave States 



T/T, 


1.1 1.3 1.5 2 2.3 


M[xc{lP)] 


3.38 - - - - 


Eb[Xc{IP)] 


« - 




9.95 10.05 10.11 10.23 10.30 


EB[Xb{lP)] 


0.35 0.25 0.19 0.07 w 


M[xi>{2P)] 


10.25 10.30 - 


EB[Xb{2P)] 


0.05 « - 



TABLE I: Summary of masses and binding energies (in [GeV]) for P— wave quarkonia in the QGP as 
extracted from the finite-temperature T— matrix determinant, Eq. ([7|. 

Next we study Q-Q scattering in a relative P— wave. In order to assess the formation of bound 
states in this channel we rely on the condition in Eq. ([7]) and determine the zeroes of det^(i?), 
which is plotted in Fig. [4]for several temperatures for both charm and bottom. We only find one c-c 
bound state at the lowest temperature (T = 1.1 Tc), aX E 3.4 GeV (just below threshold), which 
we associate with the IP charmonium, Xc- As the temperature increases, the Xc state rapidly shifts 
into the continuum. 

The P— wave h-h system exhibits two bound states at the lowest temperature, which we may 
identify with the Xb{^P) and Xbi'^P) as their energies {E = 9.95, 10.25 GeV) are close to the 
nominal values in the vacuum. The Xfe(2P) state moves beyond threshold for T ~ 1.3 Tc and the 
(IP) state for T ~ 2.3Tc. Both the mass and the binding energies, {Eb = Eth — M), of the P— wave 
states are summarized in Tab. Ufor several temperatures. 



3. Continuum Scattering 

The T— matrix approach also encompasses the continuum part of the spectrum. This is not 
easily appreciated in Figs. [2] and [3] because of the different scales of the narrow bound state signal 
and the amplitude above threshold. Note that the determinant of the transition matrix, det^, 
vanishes for some energies above threshold (cf. Figs. [2]|3]), possibly indicating resonant scattering 
in the continuum. Indeed, as can be seen in Fig. [Sj the T— matrix for c-c S— and P— wave 
scattering exhibits substantial correlations above threshold. The imaginary part shows a distorted 
resonant shape, which peaks at approximately the same energy where the real part vanishes. 
The non-perturbative effects of the Q-Q rescattering above threshold are evident as we compare 
the T— matrix to its Born approximation, V (dash-dotted line in Fig. (H). We confirm that the 
T— matrix shows the expected behavior at high energies, i.e., the real part converges to the Born 
approximation and the imaginary part tends to vanish. Finally, Fig. [6] displays the imaginary part 
of the S— wave scattering amplitude on a logarithmic scale over a wide energy range below and 
above the Q-Q threshold. 



4-. Sensitivity to IQCD Potential 



In Fig. [7] we show the c-c 5— wave scattering amplitude based on the potential fitted directly to 
the internal energy data of Ref. [25|]. We have kept rUc = 1.7 GeV for comparison with our previous 
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FIG. 5: c-c scattering amplitude in S— (up) and P— wave (below) above threshold. The Born approximation 
to the amplitude is also shown (dashed-dotted). 




10 12 14 

E. [GeV] 



16 18 



FIG. 6: Imaginary part of the S"— wave scattering amplitude including bound and scattering parts of the 
spectrum (left, charmonium; right, bottomonium) at T = 1.1 Tc. 



results. The T— matrix exhibits two bound states with a stronger binding, as to be expected from 
the potential comparison in Fig. [TJ With the same bare quark mass, the ground state is located 
at a much lower energy, E « 2.2 GeV, translating into a binding energy of about 1.2 GeV, in 
agreement with the results of Refs. 0, |20] within a Schrodinger equation using the same potential. 
As the temperature increases, the bound states rapidly shift to higher energies, reflecting the rapid 
reduction of the potential strength at low temperatures. This trend slows down beyond 2Tc, 
and the ground state eventually dissolves at T ~ 2.5 Tc. The large binding within this potential 
requires appreciable bare quark masses (ttt-c ~ 2 GeV) in order to reproduce the nominal position 
of the charmonium ground state in the vacuum. The strong attraction is presumably related to the 
entropy contribution to the Q-Q free energy at large distances (cf. Figs. 3 and 4 in Ref. |25j), which 
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peaks at Tc and decreases steeply with temperature. The large-distance limit of the internal energy, 
Ui°, which inherits this behavior, is subtracted to generate the Q-Q potential, cf. Eq. (jlOp . As 
mentioned in Sec. IIII Al might be interpreted as a contribution to the in-medium quark mass, 
m*(T) = mc + Uj^{T)/2, i.e., a quark selfenergy contribution. However, if no further r-dependence 
(or momentum dependence) is considered, the simple subtraction of from the internal energy 
distorts the normalization of the potential at short distances, where it should be described by 
perturbative QCD (one-gluon exchange). A more detailed investigation of these interplays will be 
carried out in future work. 




FIG. 7: Same as in Fig. [5] for the potential derived from the IQCD internal energy of Ref. [25]. From left to 
right and up to down the temperatures are (1.10, 1.15, 1.20, 1.50, 2.00, 2.50) Tc- 



IV. QUARKONIUM SPECTRAL FUNCTIONS AND EUCLIDEAN-TIME CORRELA- 
TORS 

A. Spectral Function 

The T— matrix formalism used above can be directly applied to evaluate mesonic spectral func- 
tions for the different quarkonium channels. The spectral functions encode the information on both 
the bound and scattering states in the continuum (E > Efh), similar to the T— matrix. Moreover, 
they allow for a quantitative connection between the present approach and Euclidean-time correla- 
tion functions, which have been calculated in lattice QCD with rather high precision a a. Such a 
comparison has recently been conducted in Ref. 0] where the heavy-quark interaction in the QGP 
has been studied by solving the bound-state problem using a Schrodinger equation with either 
a screened Cornell- type potential or IQCD-based internal energies (similar to the present work). 
The quarkonia spectral functions were then composed of J-function like bound states with weights 
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determined by the decay constant of the state and a continuum assuming free quark propagation 
with a threshold behavior taken from perturbative QCD sO, 31|. 



(12) 



The decay constants are related to the (derivative of the) radial wave function at the origin for 5— 
(P— )wave states [s^], while the functional form of the continuum threshold, given by f{u>,Etii), 
depends on the specific channel (pseudo-/scalar, axial- /vector) [s^, 31 1. The threshold energies 
were set to Ef^ = 4.5 GeV and E^j^ = 11 GeV, based on the phenomenological observation that no 
narrow mesonic resonances appear in the spectrum beyond this energies. The resulting correlation 
functions qualitatively reproduced the features observed in IQCD for the scalar channel {xc,b), 
whereas sizable discrepancies were found for the pseudoscalar and vector channels (r/c,fe and J/ip, 
T). 

In the present approach, the Q-Q system is interacting also above the Q-Q threshold with 
the same potential that generates the bound-state solutions, which, in particular, accounts for 
the transition between the discrete and the continuum part of the spectrum. Nonperturbative 
effects play an especially important role when reduced binding energies drive states toward and 
across the two-particle threshold. The mesonic spectral function is given by the imaginary part of 
the heavy-quark two-point (current-current) correlation function in momentum space, G{E,P), as 
pictorially depicted by its perturbation series in Fig. [HI The correlation function can be calculated 
from the T— matrix by closing the external legs with the appropriate momentum integrations and 
the corresponding current operator. Schematically, one has 



G = G^ + G^TG^ 



(13) 



where G^ is the lowest order correlation function, which represents the uncorrelated Q-Q propagator 
in a given mesonic channel, 

G'iE,P = 0;T) = iNfN, / j^^TrlTM A+{k)TM A-i-k)} GQQiE;k) [1 - 2 (iJk)] ,(14) 



(27r)3 



with A±{k) = {iOk^^ — k^ ± mQ)/2mQ the positive/negative energy projectors, = 
(1,75,7^^,7^^75) and Nf{Nc) the number of flavors (colors). We take A'^^ = 1, A'c = 3 as in 
to ensure the same normalization of the lowest order correlation function. Eq. (|14p denotes the 
finite-temperature result, and we have used the explicit decomposition of the single particle propa- 
gator, Sq, in terms of energy projectors, which recovers the BbS 3D-reduction scheme used in the 
calculation of the T— matrix. 
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FIG. 8: Diagrammatic representation of the Q-Q correlation function. The solid dots represent Ym operators 
specifying different mesonic channels. 
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The Q-Q rescattering is encoded in the second (two- loop) term of Eq. ([T3]). A proper connection 
has to be made between a potential description of the interaction and the relativistic invariant 
amplitude entering AG (see the Appendix for details). In particular, we have thus far suppressed 
the tensor structure of T (and V) in Dirac space. At high energies the Q-Q interaction should 
correspond to perturbative one-gluon exchange, which has a vector structure, whereas at low 
energies IQCD finds the potential to be compatible with a scalar structure [33]. In absence of further 
information, particularly for the intermediate energy regime, we consider both tensor structures 
alternatively and write for the matrix elements To = uT uT vT v, with F = 1,7^ and u{v) the 
positive (negative) energy Dirac spinors. This leads to the following traces to be evaluated in 
Eq. (US]): 

Tv{r M,r) = TT{A+{k)rM A^{-k)t A-{-k')rM A+{k')t} . (15) 

It turns out that they can be written in a partial wave expansion as performed for the T-matrix, 

Tr(rM,f) = al"'^{k,k')Po{cosekk') + a\'''^{k,k')P^{cosekk') + al'''^{k,k')P2{cosekk') , (16) 

so that all angular integrations can be done analytically by using the orthogonality of the Legendre 
polynomials. We thus have 

AGiE-T) = NfN,^ I dkk^GQQiE;k)[l-2fQiu;k)] 

X J dk'k'^GQQ{E;k')[l-2f'^{u;k')]T{TM,r;E;k,k') , (17) 

with the kernel T given by 

T{TM,r;E;k,k') = J d{cosekk')TT{TM,r;k,k' ,dkk')T{E;k,k') 

= 8TT[ao{k,k')To{E;k,k') +ai{k,k')Ti{E;k,k')] , (18) 

and ai coefficients as tabulated in Table |TT1 We note that for a given channel, e.g. pseudoscalar, 
in principle both the S— and P— wave components of the T— matrix contribute to the correlation 
function, whereas the usual spectroscopic (nonrelativistic) characterization of quarkonium states 
is based on orbital angular momentum quantum numbers {LS scheme). The (undesired) mixing of 
S— and P— wave components in the correlation function is related to the use of the JM (helicity) 
basis of the Q-Q spectrum at high energies, in which a different partial wave decomposition of 
the T— matrix follows. However, for the scalar and pseudoscalar channels the coefficient in Ta- 
ble HI] corresponding to the "natural" partial wave is leading in the non-relativistic (heavy-quark) 
expansion, whereas the other one, introducing an admixture of the "unnatural" partial wave, is of 
higher order, cf. Table lllli For simplicity, we shall work with the non-relativistic approximation 
for the ai coefficients, which filters the appropriate partial wave for the scalar and pseudoscalar 
channels (and for consistency with the spin-averaged nature of the interaction potential we shall 
consider pseudoscalar /vector, scalar /axialvector degeneracy for the correlation functions as done 
for the T— matrix)^. We comment below on the accuracy of this approximation. 



^ Note that at this level of approximation the distinction between a scalar- or vector-like structure for T is immaterial 
since both have the same heavy-quark limit for the o; coefficients (modulo global signs). 
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^ T 


PS,V 
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TABLE II: a/ coefBcients in a partial wave basis up to L = 1. 



rA/,f 


ao(fc,A;') 


ai{k, k') 


s,s 


0(fcVm2) 


-2m, + 0{e/ml) 


s,v 


o{e/mi) 


2 kk' 


PS,S 




0{k/m^) 


PS,V 


-2 + 0{k^lml) 






TABLE III: Lowest order of ai coefficients in a (1/m) expansion. 



B. Euclidean-Time Correlation Functions 



The Euclidean-time correlation function is defined as the thermal mesonic two-point correlation 



function in a mixed Euclidean-time- momentum representation [30, l3ll . |34| |. It can be expressed 
in a spectral representation as an integral transformation of the mesonic spectral function (here 
P = 0), 

G{t,T)= / dwa(w,r)/C(r,cj,r) , (19) 

JO 

where the kernel of the transformation, 

cosh[a;(r-/?/2)] 

iC{t,u;,T) = . , (20) 

smh(a;p/2) 

is symmetric with respect to r = (3/2 (and r € [0,/?]). The Euclidean-time correlation function 
scans the full spectrum of the system. In particular, for r ^ the kernel decreases rather slowly 
with energy and thus the correlation function is dominated by contributions of the Q-Q continuum. 
On the other hand, for r (3/2, the kernel exhibits the maximal decrease, so that the correlation 
function becomes mostly sensitive to the contribution from the low-energy region of the spectrum, 
in particular the bound states. To isolate the medium effects on the mesonic spectral function from 
the temperature dependence introduced by the kernel it has been proposed [3, [H] to normalize the 
correlation function at a given temperature to a so-called "reconstructed" correlation function, 

Gr{T,T)= duja{uj,T = 0)}C{t,uj,T) , (21) 
Jo 

which is obtained by replacing a{uj,T) by a reference spectral function (for instance the vacuum 
spectral function) and transformed with the same finite-temperature integral kernel. As in Ref. 
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we will first assume (j{uj, T = 0) of the form in Eq. (jl2p with the vacuum input for the bound- 
state part and a shape function f{ui,Eth) given by the perturbative QCD continuum. For the 

open-charm (bottom) threshold, we consider E^'f^'^^^^ = 2M^(^) = 3.74(10.56) GeV, but also 
check the sensitivity to changes in the pertinent free open heavy-flavor meson thresholds by using 
Ell^''''^ = 4.5(11) GeV as in Ref. [9]. As we shall see, the use of a simplified spectral function 
as in Eq. (|12p may introduce spurious features in the normalized correlation function which could 
mask the actual effect of the medium-modified Q-Q interaction. In order to enable a more direct 
comparison to IQCD evaluations, we will also normalize our results to actual spectral functions 
calculated in our approach. 

It is clear from Eq. (jl9p that the full energy regime of the spectral function figures into the calcu- 
lation of the Euclidean-time correlation function. The approximations introduced in Sec. IIV Al are 
expected to be reliable up to energies above the Q-Q threshold, where non-perturbative effects from 
the Q-Q interaction prevail, as is already manifest from non-trivial structure in the T— matrix. For 
higher energies in the continuum region we do not expect these approximations to hold. However, 
the high energy part of the continuum should be only relevant for r ^ 0, where the normalized 
correlation function approaches 1 and is not sensitive to the evolution of the quarkonia states with 
the temperature. 

C. Numerical Results 

1. Constant Heavy- Quark Mass and Small Width 

Following our studies of the T— matrices, we will first investigate mesonic spectral functions and 
normalized correlation functions for a constant quark mass, and therefore the continuum threshold 
is not dependent on the temperature. The S— wave charmonium spectral function is shown in the 
left panel of Fig. [9] for several temperatures, together with the uncorrelated (perturbative) two- 
particle continuum, Eq. (|14p. As expected, the spectral function exhibits the same charmonium 
bound states as found in the T— matrix, as well as their evolution to higher energies as the tem- 
perature is increased. At all temperatures, the (non-perturbative) rescattering of the Q-Q system 
dynamically generates a substantial enhancement of strength above the c-c threshold relative to the 
uncorrelated two-particle continuum (at T = 1.1 Tc the remnant of the first excited state {ip{2S)) 
is still visible). This important effect reflects the Q-Q correlations already found for the T— matrix 
above threshold (especially when a bound state passes into the continuum). 

The corresponding normalized Euclidean correlation function for the same set of temperatures is 
displayed in the right panel of Fig. [9l The normalized correlator converges to unity at r — > 0, which 
reconfirms the correct normalization of the continuum part of the spectrum (it is also symmetric 
with respect to [3/2). At low r, where the integral in Eq. (jl9p is mostly dominated by the continuum 
region, the normalized correlator moderately increases, reaches a maximum, and then rapidly drops 
for r approaching (3/2, indicating a loss of strength of the finite temperature correlator relative to 
the zero-temperature one in the low-energy part of the spectrum. The temperature evolution of the 
correlation function is a combined result of a decrease in binding energy of the bound states and 
the contribution of the non-perturbative continuum. The sizable drop at large r is in qualitative 
agreement with the IQCD charmonium S*— wave correlator [3]. The latter, though, exhibits a 
smaller reduction and a weaker r— dependence, showing appreciable deviations from unity only 
beyond T = 1.5 Tc. 
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FIG. 9: Left panel: Imaginary part of the correlated two-particle propagator (Green's function) for c-c 
S"— wave scattering at several temperatures and constant quark mass, TOc=1.7GeV. The imaginary part of 
the uncorrelated propagator (dashed line) is also shown for reference. Right panel: Corresponding normalized 
mesonic correlation functions. 



The charmonium spectral function and normalized correlator for P— wave scattering is displayed 
in Fig. [lOl As was already discussed in Sec. IIIIBl we only find a single bound state (xc) just below 
the threshold, which rapidly melts into the continuum as the temperature increases. Consequently, 
a sizable threshold enhancement effect is observed. Despite the fast melting of the P— wave state, 
the correlation function steeply rises in the low-r regime, due to (i) the contribution from the 
non-perturbative rescattering above threshold, and (ii) a larger threshold energy in the schematic 
vacuum spectral function entering the "reconstructed" correlator. For r — > /3/2, the correlator 
levels off well above unity, due to the absence of an energy gap between the P— wave state and the 
continuum, which renders the (enhanced) continuum contribution to the correlator dominant even 
for r — >■ /3/2. The main features of our results at a given temperature are qualitatively in line with 
the IQCD P— wave correlators [3]; however, the temperature dependence is not: our correlators 
attenuate with temperature whereas the IQCD correlator increases. This appears to be a rather 
direct indication that the in-medium c-c threshold is lowered with increasing temperature. 

While our results are qualitatively similar to those of Ref. P], the following observations are 
in order. In Ref. 0], the increase of the correlator at low and intermediate r is induced by the 
temperature-dependent decrease of the continuum threshold. Thus far, we have not considered 
the temperature effect on the Q-Q threshold. However, the non-perturbative enhancement in the 
spectral function around threshold and above, which is not included in Ref. [o^], turns out to be 
essential for a quantitative assessment of the quarkonium correlation (and spectral) functions, 
especially when the energy gap between the discrete and continuum parts of the spectrum is small 



or absent [26l . l35l . l36l |. Nevertheless, our results leave room for a lowering of the Q-Q threshold 



energy, since a downward shift of strength in the spectral function would improve (i) on the large- 
r— decrease in the S— wave correlator, and (ii) on the temperature dependence in the P— wave 
correlator, as discussed below. 



The bottomonium spectral and correlation functions follow a similar pattern as for the char- 
monium system. The results for 5"— wave scattering are displayed in Fig. [TTl At large Euclidean 
time, the correlator moderately decreases with temperature as the two excited bottomonia dis- 
sappear into the continuum. The contribution of the T(15) state, which survives up to rather 
high temperatures, makes the correlator fall more slowly than in the charmonium case. We find a 
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FIG. 10: Same as in Fig. [9] for c-c P— wave scattering. 

similar enhancement at low r, partly due to the continuum contribution with respect to the zero 
temperature case, which becomes more relevant as the bound-state binding energies decrease with 
temperature. 

100 r-: 1 1-1 rjt r-| 1 1 1 . 1 1 i 1 1 1 1 1 1 1 1 1 1 




E^„, [GeV] x [fm] 

FIG. 11: Same as in Fig. [9] for b-b S'— wave scattering. 



The P— wave bottomonium correlator, Fig. [T2l shows a large enhancement at all r, even larger 
than that of the P— wave charmonium correlator. A similar enhancement is observed in the scalar 
bottomonium correlator from IQCD [5]. As the temperature increases, the two bound states 
gradually move to higher energies and the correlator is notably attenuated. Again, as the ground 
state approaches the continuum the nonperturbative threshold strength in the spectral function 
is the decisive source of the remaining correlator enhancement. This reiterates the point that a 
comprehensive description of the Euclidean-time correlators should account for both the bound 
state evolution and Q-Q correlations above threshold. 



2. Sensitivity to the Reconstructed Correlator 

Although our results thus far roughly reproduce some of the trends of the normalized correlators 
from IQCD, the 5— wave correlator exhibits a marked r dependence which is inconsistent with IQCD 
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FIG. 12: Same as in Fig. [9] for b-b P— wave scattering. 



1.6 



~i ' r 



■ 




= 4.5 GeV 


• 




= 2M„ 


♦ 




= 3.4 GeV 



J 1 L. 



_1_ 



0.2 0.3 



1.0 

0.0 



■ 




= 4.5 GeV 


• 


E«, 


= 2M„ 


♦ 


E,„ 


= 3.4 GeV 



■ .. 



0.2 0.3 
T [fm] 



FIG. 13: Normalized correlation function for c-c S'— wave (left) and P— wave (right) scattering at T = 
l.lTc- As indicated in the legend, three different continuum threshold energies have been used in the zero 
temperature spectral function: E^j^ = 4.5 GeV as in Ref. the open-charm threshold, and E^^ — 2mc 
with nic — 1.7 GeV. 



results. In particular, in the low- and intermediate-r regime, our S*— wave correlators significantly 
increase even for temperatures just above Tc, while in IQCD they are essentially unmodified until 
about 1.5 Tc. As discussed above, one reason for this rise is the threshold mismatch between the 
finite-temperature spectral function and the vacuum one, which is modeled according to Eq. (I12p . 
We therefore investigate two alternative assumptions for the reconstructed correlator (which is 
used for normalization). 



First, we change the continuum energy threshold of the vacuum spectral function by varying the 
above values in the range 

^cc(bb) ^ 3 5 (10.3-11.0) GeV. In Fig.[l3]we compare the resulting 
charmonium correlation function to our earlier result (right panels in Figs. [9l and llOh at the lowest 
temperature, T = 1.1 Tc. Not surprisingly, for both S— and P— wave scattering the low-r rise of 
the amplitude is strongly reduced as the threshold energy is shifted to lower values. At large r, 
the 5— wave correlator changes only slightly, as expected since the discrete (low-energy) part of 
the spectrum in the reference spectral function is unchanged. However we note comparably larger 
changes of the P— wave correlator also at large r, since the small binding of the Xc makes the 
correlator sensitive to the continuum contribution in the entire r domain. 
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FIG. 14: Same as in Fig. [9] but using the in- medium correlator at T = 1.5 Tc as the "reconstructed' 
correlator to normalize the plotted ratio. 



Second, we make use of the observation that in the S— wave channel the (normalized) IQCD 
correlators are essentially unchanged compared to the reconstructed correlator up to temperatures 
of about 1.5 Tc- In order to remove ambiguities from the use of a simplified reference spectral 
function in the reconstructed correlator, we use instead our calculated in-medium spectral function 
at 1.5 Tc, which incorporates the simultaneous description of bound and scattering states. Thus, by 
construction, the correlator ratio at 1.5 Tc is identically one, consistent with IQCD. The nontrivial 
result is now the further temperature evolution, as displayed in Fig. [TH It turns out that the 
new normalized correlator shows a substantially weaker r dependence, decreasing with increasing 
temperature in a better agreement with the IQCD correlator. 

To summarize this section, we have shown that the correlator ratios are quite sensitive to 
the underlying "reconstructed" correlators used for normalization, translating into appreciable 
variations in the absolute magnitude of the normalized correlators (especially in the low-r regime). 
We thus find that current discrepancies between IQCD correlation functions and potential-model 
approaches (such as in the previous Section, as well as in earlier works [^) do not necessarily 
invalidate the latter, especially if one recalls the present additional uncertainty of how to define a 
Q-Q potential from the lattice free/internal energy [tI. [26|. 

3. In-Medium Heavy-Quark Masses and Widths 

We finally consider the effect of in-medium properties of the interacting heavy quarks. First, 
we incorporate an effective in-medium quark mass as extracted from the large-distance plateau of 
the internal energy, m*Q{T) = niQ -\- U^(T)/2, and we neglect a possible momentum dependence 
of this correction. We note that in our approach an effective reduction of the heavy-quark mass 
does not only modify the Q-Q threshold energy (Eth = 2a;q=o = 2mQ), but also figures into 
the two-particle propagator as a selfenergy contribution and therefore is iterated in the scattering 
equation. We again adjust the bare quark mass to recover the vacuum mass of the quarkonium 
ground state at 1.1 Tc. The charmonium spectral function for 5"— wave scattering (left panel of 
Fig. [T5]) shows a sizable downward shift of the bound state strength, which for higher temperatures 
is compensated by the reduction in the binding energy. As expected, the dissociation temperature 
of the ground state is smaller, T^iss ~ 2.5 Tc, as compared to the calculation with a fixed quark 
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mass. The Euclidean correlator (right panel in Fig. llSp . which we normalize to the 1.5 Tc spectral 
function as discussed in the previous section, has a rather mild temperature dependence, with an 
attenuation of the order of 10-15%, in a closer agreement with the S— wave IQCD correlator. Thus, 
for T > 1.5Tc, the T-dependent (decreasing) threshold improves the agreement with the lattice 
calculations. On the other hand, a strong T dependence as dictated by the large-distance behavior 
of the internal energy at T < 1.5Tc would not be supported the IQCD correlator (a recent analysis 
of charmonium properties in quenched IQCD shows a temperature dependent charmonium mass 
reduction of about 6% at 2Tc |37||). For more stringent conclusions the present uncertainties in the 
extraction of the heavy-quark potential, as well as the assumption of a 3-momentum independent 
mass correction (which upsets the perturbative normalization of the potentials at short distances 
and therefore represents an upper estimate of the in-medium mass effect), need to be scrutinized. 




2.8 3.0 3.2 

E [GeV] 




FIG. 15: Left: Imaginary part of the correlated two-particle propagator (Green's function) for c-c S'— wave 
scattering at several temperatures. The temperature dependent effective quark mass, m*{T), has been used 
in the calculation. Right: Corresponding normalized mesonic correlation function at several temperatures. 



We furthermore study the impact of finite quarkonium widths, which are an essential ingredi- 
ent in the phenomenology of charm and bottom in heavy-ion collisions. As suggested by recent 



analysis of heavy-quark diffusion in a QGP [20, |2l|] , as well as parton-induced break-up reactions 



of charmonia |11|], their widths are expected to be of the order of 100 MeV at temperatures around 
1.5 Tc. We can easily incorporate such effects by dressing the charm quarks with an imaginary 
selfenergy in the two-particle propagator, Eq. The results for the 5— wave charmonium cor- 
relation function in the fixed quark mass approach are shown in Fig. [TBI where a charm width 
of 50 MeV (generating charmonium widths of ~100 MeV) has been used, in comparison with the 
narrow-width limit. At 2Tc the Euclidean correlator is modified by only a few percent. Especially 
in view of other current uncertainties, the correlators are rather insensitive to phenomenologically 
relevant magnitudes of the quarkonium decay widths; this may after all not be surprising since F^ 
is very much smaller than the charmonium mass, and also appreciably smaller than the typical 
temperature. 
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FIG. 16: 5— wave charmonium correlator with a charm width of 50 MeV, as compared to the narrow quark 
calculation. 

V. CONCLUSIONS AND OUTLOOK 

In the present article we have evaluated spectral properties of heavy quark-antiquark interactions 
(charm and bottom) in the Quark-Gluon Plasma within a T— matrix approach which allows for 
a comprehensive treatment of bound and scattering states. The basic interaction was taken to 
be a two-body potential which, following earlier works, has been identified with the heavy-quark 
internal energy evaluated in thermal lattice QCD. The finite-temperature T— matrices reconfirmed 
the survival of ground-state (S— wave) bound states for temperatures well above Tc as found in 
earlier calculations, up to ~ 2.8 Tc (> 3.5 Tc) for r]c/ J/ip {i]b/'^), as well as for P— wave and excited 
bottomonia. The dissociation mechanism is characterized by a bound state passing through the 
Q-Q threshold, at which point a strong reduction and broadening of the (imaginary part of the) 
T— matrix is observed. Nevertheless, the Q-Q system remains strongly correlated in the continuum, 
indicated by resonant-like structures in the T— matrix (albeit at much reduced magnitude) which 
deviate from the Born approximation up to energies of 1-2 GeV above Q-Q threshold. 

We have proceeded to calculate Q-Q current-current correlation functions, which, in the timelike 
regime, follow from the T— matrix by folding with the Q-Q propagator. The imaginary part of 
the correlation function (spectral function) confirmed the importance of threshold effects in the 
dissolution mechanisms, resulting in large (nonperturbative) enhancements over the perturbative 
form of the continuum. These, in turn, give substantial contributions to the Euclidean correlators 
primarily at large and intermediate time, r, which cannot be neglected in quantitative comparisons 
to (and interpretations of) lattice QCD results. Assuming constant heavy-quark masses, some 
qualitative features of the IQCD correlators (large-r decrease in the S— waves, overall increase in 
the P— waves) could be reproduced. However, the magnitude of the signal in the 5"— wave channels, 
as well as the temperature dependence in the P— wave channels, were inconsistent with IQCD. 

We have found an appreciable sensitivity of the Euclidean correlator ratios to the so-called 
"reconstructed" (reference) spectral function used for normalization (usually taken as a vacuum 
form). The choice of shape and onset of the continuum introduces r— dependencies in the normal- 
ized correlator, which may render the identification of medium effects in the calculated spectral 
functions more difficult. E.g., when normalizing to a reference spectral function calculated at 1.5 Tc 
for the S— wave c-c system (as suggested by the corresponding IQCD result), the IQCD correlators 
for higher temperatures can be much better reproduced. We have also inferred that the IQCD cor- 
relators favor a temperature-dependent decrease of the heavy-quark mass {Q-Q threshold), which 
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pushes spectral strength to lower energies and improves on the large-r and temperature dependence 
of the correlators. 

Our approach furthermore allows for establishing a closer connection to quarkonium phe- 
nomenology in heavy-ion collisions by incorporating finite width effects. E.g., when implementing 
in-medium heavy-quark widths of ~50 MeV in the two-particle propagator of the scattering equa- 
tion (inducing a charmonium width of ~100 MeV as suggested by phenomenology), we find only 
few-percent changes in the Euclidean correlators, which are superseded by other current uncertain- 
ties. 

In conclusion, our results suggest that IQCD-based potential approaches, when consistently im- 
plementing both bound and scattering states in a non-perturbative scheme, are a viable means to 
quantitatively interpret the rather precise IQCD computations on Euclidean correlation functions, 
and thus evaluate the properties of quarkonium spectral functions in the QGP. Significant uncer- 
tainties still reside in the extraction of an appropriate Q-Q potential, as well as in the determination 
of the in-medium open-charm and -bottom masses. Once quantitative agreement between model 
calculations and IQCD correlators has been established, applications to high-energy heavy-ion col- 
lisions will subject the theoretical results to experimental tests. Hopefully this will pave the way to 
progress on the long-standing challenge of connecting heavy-quarkonium observables to properties 
of the finite-temperature QCD phase transition. 
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APPENDIX A: 3-DIMENSIONAL REDUCTION OF THE BETHE-SALPETER EQUA- 
TION 

The Bethe-Salpeter equation for Q-Q scattering with the Blankenbecler-Sugar (BbS) three- 



dimensional reduction of the two-particle propagator [15l . |38| reads, in the CM frame, 



M{E,q',q)=V{q\q)+ I ^V{q\ k) ^ ^ \ ^ M{E- k, q) , (Al) 

where the invariant amplitudes M. and V are actually operators (truncated amplitudes) which act 
in the direct product of the Dirac spaces of each fermion. The BbS scheme, originally formulated for 
the nucleon-nucleon {NN) interaction, exploits the following decomposition of the single-particle 
propagator in terms of positive- and negative-energy states, 

SF{k\ k) = . - . (A2) 

Consequently, the full four-dimensional two-particle propagator, iSQ{k + P/2)SQ{k — P/2), is 
replaced by the following function 

Hk")^^#^, (A3) 

which has the same discontinuity across the right-hand cut and puts the quark (anti-quark) on the 
positive (negative) energy shell, suppressing virtual anti-quark (quark) contributions. Note that 



22 



FIG. 17: CM kinematics of the Bethe-Salpeter equation. 



since both fermions are equally off-shell the energy transfer at the interaction is zero (BbS neglects 
retardation effects), and this allows for a description in terms of a static potential, V{q',q), as 
done for instance in the Boson Exchange Model of the interaction 38| . 

One can take matrix elements in Eq. (lAlh between the appropriate Dirac spinors, T[V] = 
u{q) v{—q')Ai[V]u{q')v{—q) (see Fig. [T7] for kinematics), and then Eq. ()Aip can be rewritten as 
(helicity indices omitted) 

T{E- q', q) = V{q\ q) + / j:^V{q', k) ^(^^ ^' ^ ' (A4) 

where we have used the following representation of the energy projectors 



A+(fc) 



2 m 



A^Ck) = . (A5) 

2 m 

The connection between T, V and the actual (static) potential V and the T— matrix in Eq. ([1]) 
can be derived by considering a tensor structure for V and performing a non-relativistic reduction 
of the resulting amplitude (for instance consider V given by the Yukawa scalar-meson exchange 
amplitude, which can be fully derived from the Lagrangian, Cs = 93 ^^4') ■ It turns out that 
\ Q) = Q) + 0{q^ /m?), with V related to the corresponding potential in coordinate space 

by 

V{r) = — ^ / (fke^^^'Vik) , (A6) 



(27r) 

and k = q' — q. The partial wave decomposition of the potential (and of T) is given by 

V{q',q)=47TY,{2l + l)Vi{q',q)Pi{coseg,g) , (A7) 

I 

and then Eqs. (Il|4p follow. The temperature dependence is accounted for by introducing a [1 — 2/*^] 
factor for each two-particle loop, and the quark selfenergy enters the two-particle propagator by 
the replacement 

[s/A-ujl + ie]-^ — > [s/i-ujl-2iujklni^-^ (A8) 
in Eq. ()A4p . with uj^ satisfying the dispersion relation in Eq. ([3]). 
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